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1 Introduction 

Much attention has recently been devoted to the perturbative "hard", or BFKL, pomeron 
[1], especially in relation to the study of the small x behaviour of the deep inelastic scattering 
structure functions (see a recent review in [2] ) . In application to soft phenomena, the value of 
the pomeron intercept is of principal importance. For the BFKL pomeron it is considerably 
above unity: 

oibfkl(0) = 1 - (Sa s /ir)E , 

where the "energy" Eq is equal to —4 In 2, and a s is the (fixed) QCD coupling constant [1]. 
As a result, to obtain a unitary amplitude one has to take into account more than two, in 
fact, any number of interacting reggeized gluons. This problem simplifies in the large N c 
limit, when it reduces to summing all multipomeron exchanges [3]. 

For the negative signatured amplitude the lowest order contribution comes from the ex- 
change of an odderon, a state formed from three reggeized gluons in a symmetric colour 
configuration [4]. Its relative importance is controlled by the odderon intercept. Should it 
also lie above unity, unitarization would require summing any number of exchanged odderons 
as well. 

It has not been possible to obtain a complete solution of the odderon equation ("the 
BKP equation") up to now. Certain encouraging ideas have been proposed, however, in [5-7], 
based on the conformal symmetry of the equation and the Yang-Baxter technique. In [8] 
the conformal symmetry was used to reduce the problem to an one-dimensional equation. 
Variational calculations on the basis of this approach, with a relatively simple trial function 
containng two free parameters, gave an intercept above unity, although lower than for the 
pomeron [9]: 

cwd(0) = 1 - {3a s /7r)E odd , E odd < -0.37. 

A different scheme of variational calculations of the odderon intercept was adopted in 
our paper [10]. We studied the odderon energy matrix in a non conformally invariant basis 
of functions, whose number was taken rather large (up to 3368), which corresponds to a 
correspondingly large number of variational parameters. Our best result for the energy 
however turned out positive 

E odd < 0.45 

corresponding to an intercept below unity. In view of the variational character of the calcu- 
lations this result evidently does not contradict [9] but seems to be much weaker. For that 
reason we did not give much importance to our result at the time when it was obtained, so 
that it remained unpublished. 
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However, further study of the odderon problem has given us some motivation to believe 
that this result may be closer to reality than the one obtained in [9]. The point is that 
the equation for the odderon at rest admits a simple solution corresponding to the energy 
exactly equal to zero and the intercept j = 1. This solution is based on the so-called bootstrap 
relation in the BFKL theory [11] and we call it the bootstrap solution. One can view this 
solution as a true bound state (normalizable) . It is nodeless and possesses a maximal possible 
symmetry. Therefore, from the experience gained in quantum-mechanical problems, we can 
expect it to correspond to the ground state of the odderon, that is, to the lowest energy 
possible. 

In Sec. 2 we discuss this point in more detail giving some mathematical arguments in 
favour of this result. In view of this, our calculations aquire a better status, which gives us 
a reason to present them in this paper, in Sees. 3-5. Some conclusions are drawn in Sec. 6. 
All mathematical details are relegated to Appendices. 

2 The bootstrap solution as the ground state of the odderon 

As shown in [4] the odderon wave function in the transversal momentum space ip(qi,q 2 , 93) 
satisfies a Schroedinger-like equation 

Hip = Eip. (1) 
The hamiltonian H is a sum of kinetic terms and pair interactions between the three gluons: 

H = Ti + T 2 + T 3 + U 12 + U 23 + U 31 . (2) 
Each kinetic term is given by the gluonic Regge trajectory with a minus sign in units 2>a s /i:: 

where, with an infrared regularization provided by the gluon mass m, 

r](q) = q 2 + m 2 . (4) 

The interaction terms U are integral operators in the momentum space of the three gluons 
with an integration measure 

d = d 2 g 1 d 2 q 2 d 2 q 3 5 2 (q 1 + q 2 + g 3 ~ g) , g x 
r){qi)v(Q2)v(Q3) 

where q is the (fixed) total transverse momentum of the odderon. In this section we shall 
consider only the case q = so that q\ + q 2 + q% = will be assumed. Explicitly the kernel, 
say, of U\ 2 is given by 

tMgi, 92, 93191,92. 93) = V(Q3)S 2 (93 - 93)^i2(9i,92|9l,92), ( 6 ) 
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where V± 2 is the BFKL interaction for two gluons in a vector colour state: 

(7) 



^12(91,92^1,^2) = 



V{Qi + 92 j 



77(91 - q'i) 

(Note that this interaction is twice smaller than the one for the vacuum chanel which enters 
the standard pomeron equation). 

It is well-known that like the pomeron equation, the BKP equation (1) is infrared stable, 
that is, it remains meaningful if one puts m = in (4) [4]. However, for our purpose it will 
be convenient to proceed in a different manner. Rescaling q — > q/m we eliminate m from Eq. 
(1) altogether, the resulting rj(q) having the form 

V(q) = Q 2 + l- (8) 

Thus one observes that Eq. (1) (homogeneous!) is in fact independent of m and so are the 
energy eigenvalues E, which is of no wonder, since they are dimensionless. In the following 
we shall discuss Eq. (1) in the infrared regularized form provided by 77 given by (8). It 
clearly shows that the singularities of the equation come only from the ultraviolet region 
when <7i i2 ,3 ->■ 00. 

The crucial point for our discussion is the existence of the bootstrap solution which is a 
constant: 

^(91,92,93) = -00- (9) 
As mentioned, the bootstrap identity is essential for this: 

v(q[)v(q' 2 ) 

which is a consequence of the structure of V and u in terms of 77 and valid for an arbitrary 
r](q). Applying this identity we find 



f d 2 q' 

J py^ ^ 12(91,92191,92) =u(qi)+ ^(92)- u{q!+q 2 ), (10) 



Hip B = ipo(-u(qi) - ^(92) - ^(93) 

+2u(qi) + 2u(q 2 ) + 2uj(q 3 ) - u{q 1 + q 2 ) - ^(92 + 93) - ^(g 3 + qx)) = 0. 

Here the first three terms come from the kinetic part and the rest from the interaction 
according to (10). This expression vanishes because at q = we have oj(q\ + q 2 ) = ^(93) 
etc. Evidently the solution (9) is normalizable with the measure given by (5) and 77 taken 
according to (8): 

||<M| 2 = |V'o| 2 /^<oo. (11) 

Therefore it represents a true bound state of the three gluons with a zero energy. It evidently 
possesses the maximal symmetry possible. As mentioned, our experience with quantum- 
mechanical problems then suggests that it is the ground state of the system. 
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To somewhat strengthen this proposition, we take a bit more mathematical point of view. 
The integral equation (1) is singular for two reasons. One is evident and comes from the 
pairwise nature of the interaction. It is common to all quantum mechanical three-body 
problems. As is well-known, it can be cured by going over to the equivalent Faddeev matrix 
equations for parts of the wave function. However there is another source of the singularity 
related to a bad ultraviolet behaviour of the kinetic terms: they grow very slow with q — > oo 
(only logarithmically) making the equation badly singular. This singularity persists also 
in the corresponding two-body equation for the BFKL pomeron and is responsible for its 
spectrum to be different from the "free" equation with V = 0. To cure this singularity we 
introduce an ultraviolet cutoff into rj substituting (8) by 

= v(q)eMvq 2 ), ^>o (12) 

With v > the kinetic terms grow exponentially with the momentum. One easily finds that 

2 1 

^1/(9)9-00- 2 ex P(o^ 2 )- ( 13 ) 

Now one can apply the standard methods to study the spectrum of the odderon regularized 
in the ultraviolet. 

Let us illustrate how it works with a much simpler two-body problem of the pomeron at 
rest. After changing its wave function ip{q\, 92) with q\ + c/2 = according to 

^(91,92) = Viy(qi)4>(qi) 

the pomeron equation reduces to 

(H (q) - E)<P{q) = - J d 2 q'U(q, J)tj>tf), (14) 
with Ho(q) = —2oj u (q) and the kernel U given by 

From (14) we go over to an equivalent eigenvalue equation: 

K(E) X x = \(E)xx, (16) 

where 

K(E) = (H -E)- 1 / 2 U(Ho-Er 1 / 2 . (17) 

Evidently the energy eigenvalues E for Eq. (14) are determined from the ones for (16) by 
the equation 

X(E) = -1 (18) 
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and the corresponding eigenfunctions are related by 

<t>= (Ho-Ey^x^. (19) 

Eq. (16) has usually much better properties as compared to the initial Schroedinger equation 
(14). The denominators in (17) normally make the kernel K to be of the Fredholm type, 
provided that both Hq and U are not too badly behaved and that E does not lie inside the 
spectrum of Hq. Then Eq. (16) possesses only a discrete spectrum of eigenvalues, which may 
be ordered according to their absolute values as 

I Ai| > |A 2 | > .... 

With a finite norm of K, the number of A's with an absolute value more than unity is finite. 
According to (18) this means that only a finite number of discrete eigenvalues E may exist 
below the continuum spectrum, the latter coinciding with that of Hq. 

In the pomeron case we find by a direct calculation that with v > the norm of the kernel 
K{E) is finite for all E below the minimal value of Hq: 

\\K u {E)f = J d 2 qd 2 q'\K„(E;q,q')\ 2 < oo. (20) 

One can even calculate the norm at E = and u — > oo: 

ll^(0)||^oo = 64/9. (21) 

Note that this limit corresponds to q 2 ~ 1/v and consequently to an rj(q) in the form of a 
pure exponential 

r}v{q)u->oo ^ exp(i/g 2 ). (22) 

What do these results tell us about the physically interesting pomeron spectrum at v — > 0? 
Very little indeed. One observes that with a large v there can be at most 7 eigenvalues \ 
with an absolute value greater than unity at E = 0. This means that there are at most 7 
discrete negative energy levels at large v. As v diminishes the norm of K becomes larger 
and the number of discrete negative levels E also increases. When v — > the norm blows up 
to infinity and so does the number of negative energy levels. The latter finally form a cut 
which completely fills the gap between the BFKL energy level Ebfkl = — (12a s 7r) ln2 and 
the kinetic energy threshold Eq = 3a s /2ir. Except for this behaviour with the change of u, 
no new information can be extracted from this approach in the pomeron case. 

However, applied to the odderon case this argument leads to certain important conse- 
quences. One notices that with r\ changed to r\ v the bootstrap identity (10) remains valid 
and thus the bootstrap solution persists with the energy eigenvalue E = at any v. This 
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means that for the odderon equation similar to (16) a curve Xb(E) always exists which passes 
through minus unity at E = 0: A#(0) = — 1. Since different A's cannot cross with changing v 
and since |A|'s are expected to fall monotonously with \E\ for E below the spectrum of Hq, 
we find that the curve Xb(E) divides all A's into two separate sets: those which lie above it 
and those which lie below it (Fig. 1), and that this division is conserved as v changes. Since 
we expect the norm of the operator K in the properly formulated (Faddeev) equation to be 
finite at E = 0, the number of A's below Xb has also to be finite and it is conserved with 
the changing v. Then, as a first result, we find that the solutions of the physical odderon 
equation with v = with negative energies form a discrete finite set. A second result is that 
the number of such solutions can be studied at any chosen value of in particular at large v, 
since this number is adiabatically conserved with v. The odderon theory with a large v and 
the exponential r\ of the form (22) is much simpler than for the physical odderon with v = 0. 
Direct variational estimates reveal that in this case there are no solutions with nonpositive 
energies E at large v except the bootstrap one (see Appendix 1.). This brings us to the 
conclusion that the bootstrap solution is indeed the one with the lowest energy. 

Of course, this argument is not absolutely rigorous. To make it such, one has to find 
the norm of the corresponding Faddeev kernel (a 3 x 3 skew- symmetric matrix of two-body 
scattering matrices). Before that one has to study these pair scattering matrices and show 
that their properties are no worse than those of the interactions given by (6) and (7). This 
seems realizable although rather difficult. Also, even knowing that the norm of K is finite, 
one cannot exclude in principle some bizarre behaviour of the eigenvalues A which might 
invalidate the above logic. For that reason we do not consider our derivation as a definite 
proof but rather as a strong argument in favour of the bootstrap solution to represent the 
odderon ground state. 



3 Variational calculations of the odderon ground state energy 

To perform variational estimates of the odderon ground state energy it is more convenient 
to set m = from the start, to simplify the explicit form of the Hamiltonian. To make it 
non-singular we pass to a new wave function by a substitution 

3 

y^ruv (23) 

i=l 

We shall not fix the total momentum of the odderon. Then for the new and with m = 
the metric (5) changes to 

3 

dfi = U qfd 2 qi . (24) 
i=i 
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The odderon equation (1) becomes 

3 

H^ = El[ q ^, (25) 
i=i 

The Hamiltonian can be written as [ 8 ] 

H={1/2)(H 12 + H 23 + H 31 ) (26) 

where is the BFKL Hamiltonian for gluons i and k in units 3a s /2ir, which in the limit 
m = acts on the wave function as 

3 3 
H ik i/> = J] q](\nqfql+4C)iP+ J[ q 2 [q 2 ln(r 2 k /A) q 2 + (i <- k)] ii + 2(q i + q k ) 2 ^(r ik = 0), 

(27) 

where r ik = ri — r k is the (transversal) distance between the gluons and C is the Euler 
constant. 

The solution of (25) may be found by a variational approach, searching the minimum 
value of the functional 

* = / H<Pqi1>*Hil>, (28) 

with the normalization condition 

3 

X{d 2 q^X{q^ = l. (29) 
i=i 

Using the symmetry of the wave functions in the three gluons, one can study a simpler 
functional 

r 3 

£ = (1/2) / J] d 2 q^*H 12 ^. (30) 
i=i 

The odderon energy is determined by the minimal value £3 of £ according to 

E odd = (3/2)e 3 . (31) 

Note that the minimal value e 2 of the same functional on functions ip(qi,q 2 ) for two gluons 
determine the pomeron energy: Eq = e 2 . 

We expand the wave function in a sum of products of individual gluon functions: 

3 

tp(ri,r 2 ,r 3 ) = ^ c aua2ja3 *[[if) ai (ri), (32) 

ai ,02,0:3 i=l 

where the one- gluon functions V , a( r «) form a discrete complete set and are orthonormalized 
according to (29): 

J d 2 rip* a q 2 ip a < = S a>a i. (33) 
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The coefficients c ai)Q , 2ia3 have to be symmetric in all a's by the requirement of the Bose 
symmetry and normalized according to 

l c ai, 02,03! = (^) 

"1,02,03 

The two-gluon Hamiltonian H12 acts nontrivially only on the wave functions for the gluons 
1 and 2. So the energy functional becomes 

£ ~ c oi ,02,03 c a[ ,02,03^01 ,02,0^ ,a' 2 1 (35) 



01,02,0^,02,03 



where the matrix S a ^ a ^ a i a > is the two-gluon energy in the basis formed by functions ip a . 
With this matrix known, the problem of minimization of the functional £ reduces to finding 
the minimal value of a cuadratic form, that is, the minimal eigenvalue of the matrix £ QljQ2ja ' a ^ 
considered as a matrix in independent initial and final 3-gluon states. The latter means that 
this matrix should be multiplied by unity matrix for the third gluon and then symmetrized in 
all initial and final gluons. The procedure is quite straightforward, once the basic functions 
ipa are chosen. It however involves a numerical evaluation of the energy matrix elements and 
a diagonalization of the matrix, whose dimension is rapidly growing with the number of the 
basic functions taken into account. 



4 Two-gluon energy matrix for given angular momenta 

Introducing the individual gluon angular momentum I we take 4> trivially: 

^o(r) = 4>k,l( r ) ex P ( 36 ) 

where k = 0, 1, 2, ... enumerates the radial functions. In the following, instead of r, we shall 
use the variable z = lnr 2 in most cases. In terms of z and (f> 

^VmW = "(4/r 2 )(<9 2 2 - (l/4)Z 2 )W(r). (37) 

Wave functions with different values of the angular momentum are automatically orthogonal. 
For coinciding I the normalization condition for the radial functions reduces to the standard 
form for functions 

ZkAz) = (d+\l\)i> k;l (z), (38) 

which should satisfy 

J dze k ,i(z)^,i(z) = (l/47r)<W. (39) 
We assume that the radial functions are chosen to be real. 
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With the angular dependence of the wave function explicitly given by (36), one can do 
the azymuthal integrals in the potential energy in a straightforward manner. Let Oj = 
{ki,li} and take the transition between two gluon states 0:1,012 — ► 03,04. Evidently the 
total angular momentum is conserved so that the energy matrix elements are zero unless 
h + h = h + h- According to (27) the potential energy consists of two parts, the first part 
U given by an essentially Coulomb interaction and the second one Q given by a contact 
interaction, proportional to their total momentum squared. Let us begin with the Coulomb 
part U. Its two terms evidently give the same contribution due to the symmetry under the 
interchange of gluons 1 and 2. So we can take only one of them and drop the factor 1/2. 
Denote 

Vktl (z) = (d 2 - (l/4)l 2 )^ k>l (z). (40) 
Then after doing the azymuthal integration we obtain 

U ai 16tt 2 J dz 1 dz 2 r] ai (zi)i/) a3 (zi)ip a2 (z 2 )r] a4 (z2)Ui(z 1 ,z 2 ), (41) 

where I = \h — h\ = \l 2 — h\ is the angular momentum transfer and the function Ui(z±, z 2 ) is 
given by 

^ = -(l/0exp(-(//2)|z 1 -z 2 |), Z/0, (42) 

and 

Uq = max-jzi, z 2 }. (43) 

The contact part Q involves gluonic wave functions taken at the same point. After 
performing the azymuthal integration and integrating once by parts in the variable z we 
obtain 

= 8n 2 J dz [(3 + (1/2)1/! + l 2 \)^ a2 ] [{8 + (l/2)|fe + hWcvM ■ (44) 
One can somewhat simplify this expression by noting that 

(8 + (l/2)|/i + l 2 \)lp ai 1p a2 = £ai^a 2 + V'aiCaa + ^12^ai^a 2 , (45) 

where 2A\ 2 = \h + l 2 \ — \h \ — \l 2 \ and similarly for the second factor in (44). Then finally 

= 8vr 2 J dz(^ ai ip a2 + ip ai £ a2 + A 12 ip ai il) a2 )(t a3 il) a4 + ipa 3 €a 4 + A u ip a3 ^ a4 ). (46) 

The kinetic energy is easier to calculate in the momentum space. So we transform the 
basic functions to the momentum space according to 

^a(q) = J (d 2 r/2vr)V'a(r)exp(-iqr) = (-i) 1 expilcp J rdrip k j(z)Ji(qr). (47) 
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where Ji is the Bessel function. To do the integral over r it is convenient to introduce a 
Fourier transform of the function -ip with respect to the variable z: 

A,i (z) = J (dv/y/2n)if> ktl H exp ivz. (48) 

Putting this representation in (47) and doing the r-integration we obtain 

^(q) = (2/q 2 ) exp il<j> J (dv/V^)f k ,i(v)q- 2iu , (49) 

with 

Sk,liy) = {-ifh 2w {\l\/2 + iv)$ Kl {v)T{\l\/2 + iv)/T(\l\/2 - iv). (50) 

With the gluon wave functions in the momentum space given by (49), both radial and 
azymuthal integration in q are easily done. The final matrix element of the kinetic energy T 
results as 

r Q1 , Q3 = -4ni J dvf kuh {v){d/d V )f k3:h (v), h = l 3 (51) 
The differentiation gives 

{d/du)f k3>h (v) = f k3:h (v) [2i In 2 + 2iRe ^(\h\/2 + iv) + (d/dv) In [{\l\/2 + iv)cj> Kl {v)]] . 

(52) 

Correspondingly the kinetic energy matrix element separates into terms 

iSas = 87r / ^/fc ll i 1 MAs,iiMPn2 + Ite^(|Zi|/2 + »!/)] (53) 

and 

?S«3 = -47ri J dv [(\h |/2 + iv)cf> kull (v)}*(d/dv) [(1^1/2 + ^)^3^^)]. (54) 

The function (\l\/2 + iv)<f) k j(v) is nothing but the Fourier transform of £ k ,l(z) with respect 
to z. Correspondingly we denote it as 

{\l\/2 + iv)4> Kl {v) = ^ l {v). (55) 

The part can then be written as 

Ti% 3 = -4vri J d^ kull (uy(d/du)^ 3:h (u). (56) 

The orthonormalization property (39) transforms into the analogous property in the v 
space 

J duC Kl {v)iwM = (V4tt)<W. (57) 
Noting that f k i{v)f k ij(v) = £| A v )ik' we observe that the term ln2 in (53) will add a 
constant 2 In 2 to the energy. Separating another constant term 2^(1) we finally present the 
part in the final form 

T« a3 = 2(ln2 + V>(l))<W 3 +87r / [Re ^{\h\/2 + iv) - </>(!)] . (58) 
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The first, constant, term cancels an identical one in the initial Hamiltonian (27). Using the 
representation 

roo 

ip(x) - V(l) = / dt [exp(-t) - exp(-xt)] /(l - exp(-t)) (59) 
Jo 

and the othornormalization property of the set £ a we may cast in the form 

roo 

T^ a3 = 2j Q (dt/(expt-m5 aia3 -exp(t(l-\h\/2))g aiaa (t)), (60) 

where 

9a ia3 (t) = An J di>C ai {v)t,a 3 {v)cosvt. (61) 

Note that (59) is not valid for Re x = 0. Therefore this formula cannot be applied when the 
gluon orbital momentum is zero. In this case one may use 

+ ip(—iv) = ip(l + iv) + ip(l — iv), 

which formally corresponds to changing the angular momentum to be equal to 2. 

As to the second part of the kinetic energy , it turns out to be cancelled by a similar 
contribution coming from the monopole part of the Coulomb interaction for the angular 
momentum transfer equal to zero (see Appendix 2.). 

Most of the contributions to the energy presented in this section can hardly be further 
simplified and were used in the numerical calculations as they stand. The exception is the 
monopole part of the Coulomb interaction corresponding to (41) with I = (Eq. (43)). This 
part contains contributions which cancel the term in the kinetic energy and partially the 
contact interaction contribution for I = 0. The cancellation between the monopole Coulomb 
interaction and the kinetic term is responsible for the scale invariance of the energy. 
Calculation of the monopole Coulomb part is discussed in Appendix 2. 

5 The choice of basic functions and the numerical results 

A natural orthonormal discrete basis for z = lnr 2 varying from — oo to +cx) is formed by 
the harmonic oscillator proper functions. Thus we choose functions ik,i{ z ) independent of I 
and given by 

£ k (z) = c k H k (z)exp(-z 2 /2), (62) 

where H k are the Hermite polinomals and c k are determined by the normalization condition 
(39) to be 

c 2 k = l/(4n 3/2 2 k k\). (63) 
The Fourier transformation to the v space gives 

= (-i) k c k H k {u) exp(-^ 2 /2). (64) 
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In the coordinate space the function rj k j(z) is obtained from £ by differentiation: 

m 4z) = (d-(l/2)\l\)a k (z)- (65) 
The function ip k j is obtained from £ k as a solution of the differential equation 

Uz) = {d+{l/2)\l\)^ Kl {z), (66) 
with a boundary condition ipk,l(~°°) = 0. It is given by an integral 

iPk,i(z)=[ Z dz'Uz')eM-\l\(z-z')/2). (67) 

J —CO 

With this set of functions the potential part of the energy was calculated numerically. 
As to the kinetic part, the function g entering (60) can be found analytically. For transition 
k, I — ► k',l it is equal to zero if k + k' is odd. For even k + k' = 2s and k > k' 

k' 

g kk ,{t) = ^/\-l) d c k c k , exp(-i 2 /4) £ 2*>p\C p k C^-t 2 )^ , (68) 

p=0 

where 2d = k — k'. 

After the energy matrix £ aija2>a ' a ' is calculated and properly symmetrized in the three 
gluons, its lowest eigenvalue is determined, which gives an upper limit on the exact odderon 
energy according to Eq. (31). To study the minimal energy only states with the total angular 
momentum equal to zero have been included. 

The selected set of basic one-gluon functions is characterized by the maximal value of the 
angular momentum included l max and numbers of radial functions included for each wave. As 
calculations show, best results are obtained when one raises l max and the number of radials 
in all waves simultaneously. So we present here the results for the case when the number of 
radials r is the same for all angular momenta and is equal to the number of angular momenta 
included: r = l max + 1. Such a set of functions is thus characterized by a single parameter r. 
With a growth of r the number of states N rises very rapidly. In our calculations the number 
r was limited to 6. 

To study convergence of the procedure, it was first applied to the pomeron, with only 
two gluons, where the exact energy is known. The maximal r = 6 corresponds to 201 basic 
two-gluon states in this case. For the odderon with r = 6 the number of basic states rises to 
3368. 

The results of the calculations of the ground state energies of the BFKL pomeron ei and 
the odderon e% are presented in the Table for different values of r. One observes that for the 
pomeron the obtained energies are still rather far from the exact value. Thus the convergence 
of the method is rather slow. The Table also reveals that the odderon energy is essentially 
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larger than the pomeron one for a given r. So our results confirm that, in all probability, 
the odderon intercept is lower than that of the BFKL pomeron. Our values for the odderon 
energy stay positive, approaching zero quite slowly. As discussed in Sec. 2, this seems to 
favour the minimal odderon energy being equal to zero. 

6 Conclusions 

We have presented arguments that the odderon intercept is exactly equal to unity, starting 
from the explicitly known bootstrap solution. This solution corresponds to the odderon at 
rest. However for lack of scale, the intercept should be the same for nonzero momentum 
transfers. Our variational calculations seem to confirm this result. 

If it is correct then experimentally we would expect the difference between the cross- 
sections on a particle and its antiparticle to stay constant at high energies. Also the odderon 
effects in the cross-sections themselves should decrease with energy, since the pomeron con- 
tribution rises with it. 

On the theoretical side, the intercept of the odderon being equal to one means that 
exchanges of more than one odderon are of no importance. Then the unitarization program 
essentially reduces to summing only pomerons. As mentioned this can be realized in the 
high-colour limit. However to do it for the negative signatured amplitude, one has to know 
how to couple the odderon to an external sourse perturbatively. This seems to be an unsolved 
problem. 
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8 Appendix 1. Variational estimates of the odderon ground 
state energy with an exponential tj 

With (22) the form of the kinetic part and the interaction, as well as that of the measure 
(5), simplify drastically and reduce to Gaussians. Then one can calculate the average energy 
quite easily, provided the trial function is also taken Gaussian or a sum of Gaussians. We 
performed variational estimates of the ground state energy choosing the trial function in the 
form 

4>(qi, Q2, qs) = E c i e-Mti +< & + ®, qi + q 2 + q 3 = 0, (69) 
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with variational parameters Cj and /3j > 0. 
The norm of tp turns out to be 



= C i C k N ik, (70) 
ik 

with 

N tk = (Tr 2 /3)(u + p tk r 2 , (71) 

where we have defined Pi k = + (3 k . 

The average energy in the state ip results 

(H) =J^ Ci c k E ik . (72) 

ik 

The matrix E ik is a sum of three terms: 

^ = E4?, (73) 

3=1 

corresponding to the kinetic part (j = 1), the first term in the interaction (7) (j = 2) and to 
its second term (j = 3). They are given by 

4k = (3tt 2 /8KV + PikTHZv + 3Afe) _1 , (74) 

E ik = -(3vr 2 /4)(^ + f3 ik r 1 {2u 2 + 3vf3 ik + WhT 1 , (75) 

E^ = (3n 2 /8)(u + &)"> + (3 k )-\2v + 3(3 ik y l . (76) 

Numerical calculations show that whatever the number of gaussians is taken in (69) and 
however the values of fii are chosen, the matrix E ik has no nonpositive eigenvalue unless all 
Pi = 0, which case corresponds to the bootstrap solution. Thus the bootstrap solution is the 
only one with nonpositive energy in the limit v — > oo. 

9 Appendix 2. Monopole part of the Coulomb interaction 

Explicitly the monopole term contribution is given by 

/oo rz\ 
dz 1 rj ai (z 1 )t(j a3 (z 1 )z 1 / ip a2 (z2)r] a4 (z 2 ) + (ai <-> a 4 ,a 2 <-> a 3 ). 
-oo J —oo 

(77) 

Here and in the following it is assumed that I = 0, that is, l\ = Is and l 2 = h- Introduce a 
function 

Xa 2 ,a 4 (z) = j dz'ip a2 (z')r) a4 (z'). (78) 
J— 00 

Once integrating by parts we find 

Xa 2 ,a 4 {z) = 1p a2 (z)5, a4 {z) - £, a2ia4 {z), (79) 
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where the function £a2,a 4 (£) with two indices, symmetric in these, is defined as 

£02,04(2) = [ dz'£ a2 (z')£ a4 (z'). (80) 

J — oo 

As z -»■ oo, according to (39), £a 2 ,a 4 (z) -»■ (l/47r)5 a2ja4 , so that 

Xo 2 ,o 4 (°o) = -(l/47T)5 Q2ja4 . 

Having this in mind we can rewrite (77) in the form 

/oo 
(z) + («i <-> a 4 , a 2 <->■ 03)- (81) 
-00 

Integrating by parts, the integral transforms into 

/oo 
dz{Xaa,ai(z) + (V 47 >")£oi ,03)^X02,04^) + Xo 2 ,o 4 )- (82) 
-00 

The term coming from the product Xo 3 , 01 -2X02,04 cancels the contribution (a\ <-> 04, 02 CK3) 
in (81) so that the monopole contribution becomes 

/oo 
Q?2 ,04 (.ZXce2 ,04 (^) ~^Xa2 ,04 

(z))). (83) 
-00 

Now we substitute the functions x i n (83) by the symmetric functions £ using relation 
(79). Take the the first term in the integrand of (83). With (79) we obtain for it 

X03, 01X02, 04 ^03 £01 £04 Y'o3£oi£o2,04 £0:3,01 ^02^04 "I - £o3,oi£o2,04- 

Having in mind the subsequent symmetrization with respect to the interchange of gluons 1 
and 2, we can change a\ <-> 02 and 03 <-> 04 in the second term. Summed with the third 
term it then gives 

- £03,01 (^£04 + £02^04)- (84) 

Recall now that £ a2 = (d + (l/2)|Z 2 |)Y> a2 an d similarly for £ a4 . Integration by parts allows to 
substitute (84) by 

(£o 3 £oi - \h I £03 ,01)^02^04- (85) 

So finally the first term in (83) leads to the following three contributions to the monopole 
Coulomb energy: 

~ (1) 2 f°° 

U ai , 02:03, o 4 = — 16vr / ^£o 3 ,qi£o2,o 4 ) (86) 



/oo 
^03,01^02^04 (87) 
-OO 

and 

/•oo 

^01,02 ;o 3 , 04 = - 1671 " 2 / ^(^03 £01^02 £04 +£03 £01 V'oa^)- (88) 
— 00 
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Of these terms the first is divergent in its present form. It will receive its meaning after 
adding new contributions coming from the rest of the terms in (83). For that reason we have 
denoted it with a tilda. 

Now for the rest of the terms in (83). Changing the function x by £ according to (79) we 
have under the integral 

£02,04 ~\~ ^£02,04 £02,04 

+ ztp a2 (d - (l/2)|/ 2 |)£a 4 - 

Integration by parts transforms it into 

— £02,04 _ ^£o2£o4 - (89) 

The first term can be combined with (86) to give the final U^: 

/oo 
^£o2,o 4 ((l/4vr)(5 aiiQ 3 -£a 3 ,ai)- (90) 
-00 

Now the integral converges due to the property (39). Putting here the explicit form of the 
functions £o;,o fe and integrating over z we obtain the term in its definitive form: 

f / o 1 ) ,02;o 3 ,04 = 167f2 / dz l dz 2( Z l ~ Z 2 )9{Z 1 - 22)£o 1 (2l)£o 3 (2l)£o2(^2)£o 4 (^4). (91) 

The second term in (89) gives the last contribution to the monopole energy 

roc 

koi,02;o 3 ,04 = 47T(5 Ql)Q , 3 / d,ZZ^ a2 £ a4 . (92) 

J — oo 

This term cancels with the contribution to the kinetic energy. Indeed after the Fourier 
transformation to the v space according to (48), the factor z goes into idjdv. One can then see 
that (92) gives exactly the contribution T^ 2 \ Eq. (56), with an opposite sign and with gluons 
1 and 2 interchanged, which is of no importance because of the subsequent symmetrization. 

The term cancels with the part of the contact interaction Q, Eq. (46), which does 
not contain factors A: 

roo 

,02:03,04 8 71 " / ^(^03 £01^02 £04 ~\~ £a 3 £ai ^02^04 4" V'oi £02^04 £03^01 £02 ^04)- 
J — 00 

(93) 

Summed with this part gives 

/oo 
dz(lpai£a2 ~ £01^02) (£03 ^04 + Y>a 3 £a 4 )- (94) 
-00 

This expression is antisymmetric under the interchange of the gluons 1 and 2 and does not 
give any contribution to the energy. 

So finally the only contributions which remain in the interaction for zero angular momen- 
tum transfer are U^ l \ and the part of the contact interaction (46) which contains 
factors A. 
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Table 

Calculated values of the ground state energy per gluon multiplied by 2 (Eq. 
(31)) for the pomeron (62) and odderon (e^) with different numbers r of radial 
functions and angular momenta included. 



r 


£2 


£3 


1 


0.968 


0.968 


2 


0.022 


0.605 


3 


-0.475 


0.454 


4 


-0.743 


0.379 


5 


-0.912 


0.331 


6 


-1.032 


0.298 








00 


-2.773 
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